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correspond to attractive and which to repulsive forces. 
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1. Introduction 



Lattice computations of the potential of a pair of static-light mesons (in the following also 
referred to as B mesons) are of interest, because they constitute first principles determinations of a 
hadronic force. Until now interactions between static-light mesons have exclusively been studied 
in the quenched approximation [JXJ, [2|]. Here I report on the status of an investigation with two 
flavors of dynamical Wilson twisted mass quarks. Forces are not only studied between the lightest 
static-light mesons (denoted by S), but also first excitations are taken into account (denoted by P_). 
Note that there is another ongoing study of static-light meson interactions with dynamical quarks, 
which has also been reported during this conference [|0. 



2. Trial states and quantum numbers 

2.1 Static-light mesons 

Here I consider static-light mesons, which are made from a static antiquark Q and a light quark 
Y £ {u, d}. Consequently, isospin 1=1/2 and I z € {—1/2, +1/2}. Since there are no interac- 
tions involving the static quark spin, it is appropriate to classify static-light mesons by the angular 
momentum of their light degrees of freedom j. I do not consider non-trivial gluonic excitations, 
hence j = 1/2 and j z = {—1/2, +1/2}, which is the spin of the light u/d quark. Parity is also a 
quantum number, £ {+ , — }. 

The lightest static-light meson has quantum numbers = (l/2)~ (denoted by S). The first 
excitation, which is w 400 Me V heavier, has quantum numbers = (l/2) + (denoted by PJ). 
Examples of corresponding static-light meson trial states are Qy=, and Qyj\\f\Q) for S mesons 
and Q\j/\Q.) and Qyjy$\\f\£i) for P_ mesons respectively. 

For a more detailed discussion of static-light mesons I refer to [§§]. 



2.2 BB systems 

The aim of this work is to determine the potential of a pair of B mesons as a function of their 
separation R (without loss of generality I choose the axis of separation to be the z axis). To this end 
one has to compute the energy of eigenstates of the Hamiltonian containing two static antiquarks 
<2(ri) and Q{r 2 ), ri = (0,0, -R/2) and r 2 = (0,0, +R/2), which define the positions of the two B 
mesons, and which will be surrounded by light quarks and gluons. 

These BB states are characterized by several quantum numbers. Since there are two light 
u/d valence quarks, isospin / S {0, 1} and I z € {—1 , 0, +1}. Due to the separation of the static 
antiquarks along the z axis, rotational symmetry is restricted to rotations around this axis. Con- 
sequently, states can be classified by the z component of total angular momentum. However, as 



already mentioned in section [2JJ there are no interactions involving the static quark spin. There- 
fore, it is appropriate to label BB states by the z component of the angular momentum of the light 
degrees of freedom j z £ {— 1,0, +1}. Parity is also a symmetry and, therefore, a quantum num- 
ber, & G { + , — }. For states with j z = there is an additional symmetry, reflection along an axis 
perpendicular to the axis of separation (without loss of generality I choose the x axis). The cor- 
responding quantum number is & x 6 {+, — }. When using |j z | instead of j z , g? x is a quantum 



2 



Forces between static-light mesons 



Marc Wagner 



number for all states. To summarize, BB states can be characterized by the following five quantum 
numbers: (I,I Z , \j z \,&, 3^ x ). 
I use BB trial states 

(^T)AB{Qc(n)W^\ri)) (<2c(r 2 )v4 2) ( r 2))|^), (2.1) 

where the lower indices A, B and C denote spinor indices, ^ = 70/2 is the charge conjugation 
matrix and T is a combination of 7 matrices. Note that it is essential to couple the light degrees of 
freedom of both mesons in spinor space, because these degrees of freedom determine the quantum 
number \j z \. Proceeding in a naive way by coupling light and static degrees of freedom in both B 
mesons separately will not result in a well defined angular momentum \j z \ and, therefore, will mix 
different sectors. To obtain 1 = 0, the flavors of the light quarks have to be chosen according to 
^(1)^(2) _ M j _ j M> while for I = 1 three possibilities exist, yWip( 2 ) e {uu, dd , ud + du}. BB 
trial states are collected in Table [I] together with their quantum numbers. 
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Table 1: quantum numbers of BB trial states; due to explicit isospin breaking, (/ = l,I z = 0) and 
(7= l,I z = ±1) states are not degenerate in twisted mass lattice QCD (cf. section 3) and, therefore, listed 
separately; "result" characterizes the shapes of the numerically computed BB potentials (A: attractive poten- 
tial; R: repulsive potential; SS: lower asymptotic value 2m(S); SP: higher asymptotic value m(S) +m{P )\ 
cf. section 4). 



3. Lattice setup 

I use 24 3 x 48 gauge field configurations generated by the European Twisted Mass Collabora- 
tion (ETMC). The fermion action is Nf = 2 Wilson twisted mass, 

Sf\X,Z,U] = a A Y^X(x)(D^ + i^y^)x{x) (3.1) 
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[7|, where Dw is the standard Wilson Dirac operator and % = (j[ ,X ) * s the light quark 
doublet in the so-called twisted basis. In the continuum the twisted basis is related to the physical 
basis by the twist rotation y = g 1 ^ / 2 ^ where (0 is the twist angle, co has been tuned to maximal 
twist, i.e. ft) = Tt/2, where static-light mass differences are automatically G{a) improved. The 
gauge action is tree-level Symanzik improved [§]. I use j8 = 3.9 and /i q = 0.0040 corresponding to 
a lattice spacing a = 0.079(3) fm and a pion mass nips = 340(13) MeV [g]. For details regarding 
these gauge field configurations I refer to [|H||ll"|]. 

In twisted mass lattice QCD at finite lattice spacing SU(2) isospin is explicitely broken to U(l), 
i.e. I z is still a quantum number, but I is not. Moreover, parity & has to be replaced by twisted mass 
parity ^( tm ), which is parity combined with light flavor exchange. The consequence is that twisted 
mass BB sectors are either labeled by (I z ,\j z \,&> (tm) ^i' m) ) for I z = ± 1 or by (I z ,\j z \,^ tm \ ) 
for / £ = 0. A comparison with the set of quantum numbers discussed in section 22, shows that in 
the twisted mass formalism there are only half as many BB sectors as in QCD, i.e. QCD BB sectors 
are pairwise combined. Nevertheless, it is possible to unambiguously interpret states obtained 
from twisted mass correlation functions in terms of QCD quantum numbers. The method has 
successfully been applied in the context of static-light mesons [12] and is explained in detail for 
kaons and D mesons in [13]. For a detailed discussion of twisted mass symmetries in the context 
of BB systems I refer to an upcoming publication [14]. 

When computing correlation functions, I use several techniques to improve the signal quality 
including operator optimization by means of APE and Gaussian smearing and stochastic propaga- 
tors combined with timeslice dilution. These techniques are very similar to those used in a recent 
study of the static-light meson spectrum 0, |5|] and will also be explained in detail in [14]. 

In contrast to spectrum calculations for static-light mesons [Q, ||] and static-light baryons [JTJ], 
where we have always used the HYP2 static action, I perform computations both with the HYP2 
static action and with unsmeared links representing the world lines of the static antiquarks. In par- 
ticular for small QQ separations R<2a ultraviolet fluctuations are important, which are, however, 
filtered out, when using HYP smeared links. The effect of HYP smearing is shown in Figure [[]. 
For all results presented in the following potential values corresponding to R < 2a have been com- 
puted by means of unsmeared links, while for larger separations HYP smearing has been applied 
to improve the signal-to-noise ratio. 
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Figure 1: the BB potential corresponding to \jf^\j/^ 
the HYP2 static action. 



uu, r = J} computed with unsmeared links and with 
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4. Numerical results 

The BB potentials presented and discussed in the following have been obtained by fitting con- 
stants to effective mass plateaus obtained from temporal correlation functions of trial states (|Q|). 
In twisted mass lattice QCD there are 24 independent I z = trial states (i.e. trial states not related 
by symmetries) and 12 independent /, = ±1 trial states, i.e. 36 resulting potentials, which are not 
related by symmetries (cf. Table ||). Some of these potentials are quite similar, while others are 
not. In total there are four significantly different types of potentials: two of them are attractive, the 
other two are repulsive; two have have asymptotic values for large separations R, which are larger 
by around 400 MeV compared to the other two (cf. the "result" columns of Table [p). For each of 
the four types an example is plotted in Figure ||. 
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Figure 2: examples of BB potentials as functions of the separation R. 
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To understand the asymptotic behavior, it is convenient to express the BB creation operators 



appearing in ( |2.lD in terms of static-light meson creation operators. For the potentials shown in 
Figure || one finds after some linear algebra 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



(^l)A B (<2c(ri)MA(ri)) (<2c(r2W(r 2 )) = 

= -5 t (r 1 )P_ | (r 2 )+5 4 (r 1 )P_ t (r 2 )-P_ t (r 1 )5^(r 2 )+P_ ; (r 1 )5 t (r 2 ) 

(tfttXtfl(fic(ri)BA(ri)) (Gc(r 2 W(r 2 )) = 

= -5 t (r 1 )P_ ; (r 2 )+5^(r 1 )P_ t (r 2 )+P_ t (r 1 )5 ; (r 2 )-P^(r 1 )5 t (r 2 ) 

(^75)A B (<2c(ri)«A(ri)) (fic(r 2 )Mj,(r 2 )) = 

= -5 t (r 1 )5 4 (r 2 )+5 4 (r 1 )5 t (r 2 )-P_ t (r 1 )P_ i (r 2 )+P_ ; (r 1 )P_ t (r 2 ) 

C^)A B ((2c(ri>A(ri)) (<2c(r 2 W(r2)) = 

= -j5t(r 1 )5;(r 2 )-j5;(r 1 )5 t (r 2 ) + jP_ t (r 1 )P_ ; (r 2 ) + jP_ i (r 1 )P_ t (r 2 ). 
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At large separations R the BB potentials are expected to approach the sum of the masses of the 
two individual B mesons. When considering ( JOI ) to (fO|) and Figure one can see that the two 
potentials with the lower asymptotic value (yA^y/ 2 ) = uu, Y = 75 and y^y^ = uu, T = 73) 
contain SS combinations. These are significantly lighter than the also present P P combinations 
and should, therefore, dominate the correlation functions and effective masses at large temporal 
separations. The asymptotic value of the corresponding potentials should be around 2m(S), which 
is the case. In contrast to that the other two potentials with the higher asymptotic value 
(I// 1 ))// 2 ) = uu, r = 1 and y/Wy^ 2 ) = uu, T = jo) exclusively contain SP combinations. Their 
asymptotic value is expected at around m(S) + m(.P_), which is also reflected by Figure ^. 

This expansion of BB creation operators in terms of static-light meson creation operators also 
provides an explanation, why potentials computed with different operators, but which have iden- 
tical quantum numbers, are of different type. An example is given by l/Z^y/ 2 ' = uu, F = 7? and 
^(i)^( 2 ) = uu, r = 1, both having quantum numbers (7 = l,I z = +1, |y' z | = 0, & = — , £? x = —). 
The r = 73 potential is attractive with an asymptotic value at around 2m(S), while the T = 1 
potential is repulsive with an asymptotic value at around m(S) +m(P-). From (11) and ( JO] ) 
one can read off that the static-light meson content is essentially "orthogonal": the T = 73 oper- 
ator contains SS and P P combinations, whereas the T = 1 operator is exclusively made from 
SP combinations. While the corresponding T = 73 correlator yields the ground state in the 
(/ = l,I z = +1, \j z \ = 0, = — , & x = —) sector, which closely resembles a pair of S mesons, 
the r = 1 operator mainly excites the first excitation, which is similar to an SP combination. The 
generated ground state overlap is, therefore, rather small and, consequently, very large temporal 
separations would be needed to extract the ground state potential. Presumably, the potential corre- 
sponding to the r = 1 operator has a small ground state conuibution, which contaminates the first 
excited state potential. This is supported by the observation that the asymptotic value of the T = 1 
potential is slightly lower than m(S) + For a clean extraction of this first excited state an 

analysis of a 2 x 2 correlation matrix is needed. 

From the 36 independent potentials one can also deduce a rule stating, whether a BB potential 
is attractive or repulsive. The rule is quite simple. 

A BB potential is attractive, if the trial state is symmetric under meson exchange, repulsive, 
if the trial state is antisymmetric under meson exchange. 

Here meson exchange means exchange of flavor, spin and parity. One can easily verify this rule for 
the examples discussed above: the operators (4.2) and (4.4) are symmetric under meson exchange 
and give rise to attractive potentials, while the operators (4.1) and (4.3) are antisymmetric under 
meson exchange and yield repulsive potentials. This more general rule is in agreement to what has 
been observed in quenched BB computations for SS potentials [jl|, ||]. 



5. Conclusions 



I have presented results of an ongoing computation of BB potentials. Various channels char- 
acterized by the quantum numbers (I,I Z , \j z \, & , & x ) have been investigated. The computations 
have been performed with dynamical, rather light quark masses (raps ~ 340 MeV). The results 
have been interpreted in terms of individual S and P mesons. A simple rule has been established 
stating, whether a BB potential is attractive or repulsive. 
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The statistical accuracy of the correlation functions needs to be improved. BB systems are 
rather heavy and, hence, effective masses are quickly lost in noise. At the present level of statistics 
slight contamination from excited states cannot be excluded. To this end contractions are ongoing. 

Future plans include studying the light quark mass dependence, the continuum limit and finite 
volume effects. Moreover, also BB S and B S B S potentials could be computed. To treat the s quark as 
a fully dynamical quark, such computations should be performed on Nf = 2 + 1 + 1 flavor gauge 
field configurations currently produced by ETMC Jl6|]. It would also be interesting to supplement 
the lattice computation by a perturbative calculation of BB potentials at small separations R < 2. 
Finally, one could use the obtained BB potentials as input for phenomenological considerations to 
answer e.g. the question, whether two B mesons are able to form a bound state. 
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